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Recent detections of GeV photons in a few GRBs by Fermi-LAT have led to strong constraints on the bulk 
Lorcntz factor in GRB outflows. To avoid a large 77 optical depth, minimum values of the Lorentz factor have 
been estimated to be as high as 800-1200 in some bursts. Here we present a detailed calculation of the 77 optical 
depth taking into account both the geometry and the dynamics of the jet. In the framework of the internal 
shock model, we compute lightcurves in different energy bands and the corresponding spectrum and we show 
how the limits on the Lorentz factor can be significantly lowered compared to previous estimates. 
Our detailed model of the propagation of high energy photons in GRB outflows is also appropriate to study 
many other consequences of 77 annihilation in GRBs: (i) the 77 cutoff transition in a time-integrated spectrum 
is expected to be closer to a power-law steepening of the spectrum than to a sharp exponential decay; (ii) the 
temporal evolution of the 77 opacity during a burst favors a delay between the MeV and GeV light curves; 
(iii) for complex GRBs, the 77 opacity suppresses the shortest time-scale features in high energy light curves 
(above 100 MeV). Finally we also consider GRB scenarii where MeV and GeV photons are not produced at the 
same location, showing that the 77 opacity could be further lowered, reducing even more the constraint on the 
minimum Lorentz factor. 
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1. Introduction 



The compactness problem. The short time 
scales observed in GRBs (down to a few ms) can 
be used to deduce an upper limit on the size of the 
emitting region producing 7-rays. This information 
combined with the huge isotropic 7-ray luminosities 
deduced from the measured redshifts imply huge 
photon densities. Then the simplest assumption of an 
emission produced by a plasma radiating isotropically 
with no macroscopic motion predicts that 7-ray 
photons should not escape due to 77 annihilation 
77 — > e + e~ . This is in contradiction with the 
observed GRB spectra which are non-thermal and 
extend well above the rest-mass electron energy 
m e c 2 fa 511 kcV. Observation and theory can be 
reconciled by assuming that the emitting material 
is moving at ultra- relativistic velocities This 
is mainly due to the relativistic beaming. First it 
implies that the observer will see only a small fraction 
of the emitting region: the constraint on the size of 
the source is now less severe. Second, the collimation 
of photons in the same direction reduce the number 
of potential interactions. Finally the typical 77 
interaction angle becoming small the photon energy 
threshold for pair production becomes higher. This 
theoretical context combined with the observational 
data gives the possibility to estimate a minimum 
Lorentz factor r m i n for the emitting outflow in 
GRBs [2] or directly a Lorentz factor estimate if the 
77 cutoff is clearly identified in the spectrum (see Q). 

Severe constraints on the Lorentz factor from 
Fermi-LAT observations. Since the launch of 
Fermi in June 2008, the LAT instrument has detected 
high energy photons above 10 GeV in a few GRBs. 
The observed 7-ray spectrum often remains consis- 



tent with a Band function covering the GBM and 
LAT spectral ranges without any evidence of a high 
energy cutoff that could be identified as a signature 
of 77 — > e + e~. This extension by Fermi of the ob- 
served spectral range upper bound from a few MeV 
(e.g. BATSE) to 10 GeV implies constraints on r mul 
which are much more severe than the ones obtained 
previously. In a few cases r m j n has been estimated to 
be of the order of 1000 (for example: GRB 080916C - 
r min = 887 0, GRB 090510- r min = 1200 @). These 
extreme values put severe constraints on the physics 
of the central engine which should be able to strongly 
limit the baryon load in the outflow. 
However these r m i n values were obtained from a sim- 
plified "single zone" model where the space and time 
dependencies are averaged out. The motivation of this 
work is to develop a detailed approach taking into ac- 
count a more realistic treatment of the dynamics. 



2. Computing the 77 optical depth 

General 77 opacity formula. The 77 opacity (r 77 ) 
is given by: 



dl / dfl 



T 11 {E GeV ) 

l e J J E a (E GcV ,-4,) 

dE nn(E)a^(E,ip)(l - costp) (1) 

All the physical quantities are measured in the lab- 
oratory (or source) frame. EceV is the energy of 
the photon for which r 77 is calculated whereas E is 
the energy of the interacting field photon, ip rep- 
resents the interaction angle between the GeV pho- 
ton and the interacting photon and <r 77 is the 77 
interaction cross-section between these two photons. 
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Figure 1: Opacity in a single pulse: comparison 
with the semi-analytical work of [6|. 77 absorbed 
lightcurves at four different energies e = EnE/m e c 2 are 
plotted as a function of the observer time for a single 
pulse using the prescriptions and model parameters 
corresponding to the middle panel of Fig. 9 in [f| . Our 
numerical calculation (solid line) is compared to the 
semi-analytical result of [y| (dashed line). Notations are 
the same (the observer time T and observed fluxes 
eF e /Fo have normalized values). The agreement is 
excellent except for T < 10~ 2 , where the discrepancy is 
due to numerical resolution limitations (this corresponds 
to a true observer time t b s < 0.1 ms). [figure from [3]] 



E c = 2(m e c 2 ) 2 /[E GeV (l — cosip)] is the energy thresh- 
old of the field photon above which 77 annihilation 
can happen. Finally uq is the photon field distribu- 
tion [ph • cm~ 3 • erg -1 • sr^" 1 ] at a given location and 
time. 

The equation (fTJ) is made of a triple integral : the 
e?Z-integration is done over the path of the GeV 
photon from its emission location to the observer, the 
<if2-integration is done over the solid angle distribu- 
tion of the interacting photon field surrounding the 
GeV photon whereas the disintegration is done over 
its energy distribution. The equation ([I} is general 
and can be applied to any photon emitted at a given 
location and time with a given propagation direction 
within the GRB outflow. 

Validation of the model. The kernel of our study is 
the calculation of the 77 opacity created by a spherical 
flash, i.e. an instantaneous flash of photons emitted 
by an expanding relativistic spherical front. It is then 
possible to model the case of a propagating radiat- 
ing spherical front (representing for example a shock 
wave) by the succession of many spherical flashes. One 
of the critical step is the exact calculation of the pho- 



ton density nn taking into account all the relativistic 
effects. Before dealing with more complex dynamical 
configurations within the internal shock framework, 
the validity of our numerical approach was tested on 
a simple single-pulse case with a comparison to the 
previous semi-analytic study of @ (see Fig.Q]). 

3. Application to Internal Shocks 

3.1. Internal shocks within a relativistic 
outflow 

Now the model is applied to dynamical evolutions 
expected in the internal shock framework, where the 
whole prompt 7-ray emission is produced by electrons 
accelerated by shock waves propagating within a rela- 
tivistic variable outflow. We model the dynamics via 
a multiple shell model where the successive collisions 
between shells mimic the propagation of shock waves 
Each collision produces an elementary spherical 
flash: the simulated light curves are the result of the 
sum of all flashes. For each high energy photon, the 
77 opacity is computed by integrating equation (1) 
from its emission location to the observer taking into 
account the exact radiation field nn produced by all 
the collisions in the outflow. A previous study of the 
77 opacity in internal shock was made by [9(. How- 
ever the prescription used to compute r 7T was still ap- 
proximate, using the local physical conditions of the 
outflow where the high energy photon is emitted and 
applying them to an average formula of r 77 (as can be 
found in 0,111). 

3.2. Minimum Lorentz factor in GRB 
outflows - The case of GRB 08091 6C 

The first natural application of our model is 
the estimate of the minimum bulk Lorentz factor 
Turin in GRB outflows, obtained from the constraint 
T 77 (£ ; HE,ma X ) ^ 1, where EHE.max is the highest pho- 
ton energy detected in the burst. To illustrate this as- 
pect with an example, we applied our approach to the 
case of one of the four brightest GRBs detected in the 
GeV range by Fermi, i.e. GRB 080916C. The results 
are shown in Fig. [21 Using our numerical model, a syn- 
thetic GRB was generated, which reproduces the main 
observational features: the total radiated isotropic 7- 
ray energy {E iso = 8.8 x 10 54 ergs between 10 keV 
and 10 GeV), the spectral properties (E p , a, (3 pa- 
rameters of the Band function , the envelop of the 
light curve and a short time-scale variability of 0.5 s in 
the observer frame. The study is focused on the most 
constraining time bin (time bin 'b'), during which the 
highest observed photon energy was i?HE,max = 3 GeV 
(16 GeV in the source rest frame): for this reason, 
only time bins 'a' and 'b' are reproduced in the syn- 
thetic GRB. These two intervals correspond to 32 % 
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Figure 2: Minimum Lorentz factor for GRB 080916C. The two first panels are plotted for the limit case leading 
to t 77 (3 GeV) = 1 in time bin 'b', i.e. for a mean Lorentz factor F = r m i n = 340. Upper left panel: initial Lorentz factor 
distribution in the outflow. Upper right panel: 7-ray lightcurves in the GBM band (8 keV - 5 MeV, top) and at 3 GeV 
(bottom). The lightcurves are plotted as a function of t b s — iobs,tri g , where iobs.trig is the observer time of the first 
detected photons. Lower left panel: evolution of r 77 at I?he = 3 GeV against the mean Lorentz factor in the outflow L, 
following our detailed modeling (solid line) and using the average formula from 4] (dashed line). Lower right panel: 
time integrated spectrum over time bin 'b' for different mean Lorentz factors (the relative shape of the initial Lorentz 
factor distribution is kept the same) and reference spectrum without 77 annihilation (dashed line), [figure from [7|] 



of the total radiated isotropic equivalent energy. The 
minimum mean Lorentz factor r m ; n is obtained by 
requiring that r 77 (£HE,max) < 1 (see Fig. [2j lower 
panel). With the detailed calculation, we find a mini- 
mum mean Lorentz factor r m ; n = 340, i.e. a factor 2.6 
lower than the value r m ; n = 887, which was obtained 
from an approximate "single zone" model Even 
more remarkable, the whole initial distribution of the 
Lorentz factor used in this model of GRB 080916C 
(from 170 to 700) remains below the "minimum" value 



of the Lorentz factor derived from single zone models 
(see Fig. [2j upper left panel). 

3.3. Is the delayed onset of the GeV 
emission a signature of the 77 opacity ? 

The high energy emission (above 100 MeV) detected 
by Fermi in a few bright GRBs often shows a delayed 
onset compared to the softer 7-ray emission (below 
5 MeV). The analysis by [Io[ indicates that such a 
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delayed onset is present in at least 7 in a sample of 
17 GRBs detected by Fermi-L AT . This feature seems 
to be common to long and short GRB classes and its 
origin is debated 1 1 l |j . Among the proposed expla- 
nations (see e.g. (l2l Il3l. Il4l|). the possibility that this 
delayed onset is induced by a 77 opacity temporal evo- 
lution effect has already been discussed by [4J: as the 
shock wave producing the 7-ray emission expands to 
larger radii, the opacity seen by the high energy pho- 
tons evolve from an optically thick to an optically thin 
regime. The model developed in the present study is 
well appropriate to investigate this possibility in more 
details. The synthetic burst used in Fig. [2] to model 
bins 'a' and 'b' of GRB 080916C gives an example 
of a delayed onset at 3 GeV induced by an evolving 
77 opacity. The first pulse is produced at lower radii 
and in lower Lorentz factor material and is therefore 
strongly absorbed. For this reason, it is almost sup- 
pressed in the 3 GeV lightcurve, whereas the second 
pulse is well visible. Note that the model reproduces 
simultaneously the onset delay of ~ 5 s at high energy, 
and the short timescale variability of ~ 0.5 s at low 
energy. 

4. Consequences of distinct emission 
regions for MeV and GeV photons 

4.1. Are GeV and MeV photons produced 
in the same place ? 

It has been proposed in several recent studies that 
the delayed onset and/or the long-lasting tail of the 
high energy emission could be an evidence in favor of 
two different regions for the emission of MeV and GeV 
photons. An extreme version is the scenario proposed 
by [HIIIII where the whole GeV emission (prompt and 
long lasting) is produced by the external shock dur- 
ing the early deceleration of the relativistic outflow. 
Note that this scenario leads to strong constraints on 
the density and magnetization of the external medium 
[TtT ] and that the observed temporal slope of the long- 
lasting high-energy emission would imply a strongly 
pair-enriched medium (l6| . 

Even in scenarios where the prompt GeV emission 
has an internal origin, a partially distinct emission re- 
gion could be due to a spectral evolution of the prompt 
mechanism. For instance, in the framework of internal 
shocks, the evolution of the physical conditions in the 
shocked region during the propagation of a shock wave 
leads to an evolving efficiency of the IC scatterings, 
depending on the importance of Klein- Nishina correc- 
tions. This naturally leads to a variable high-energy 
component following with a delay the main (Band) 
component in the MeV range [H, [H, H3] ■ Successive 
generations of collisions in a variable outflow can also 
lead naturally to different emission regions 13]. An 



evolution in the microphysics of the acceleration pro- 
cess could also be responsible for some spectral evolu- 
tion in scenarios where there is a dominant hadronic 
component at high energy (see e.g. (HI). Finally, 
two emission regions are naturally expected in photo- 
spheric models, as it is often assumed that the main 
(Band) component has a photospheric origin and that 
internal shocks or magnetic dissipation occurring at 
larger distance produce an additional component at 
high energy (see e.g. [H, 

4.2. Loosening the constraint on r min 

As discussed in [U [25| , the possibility for the GeV 
photons to be produced in a different region than 
the MeV photons can loosen the constraint on the 
minimum Lorentz factor in GRB outflows. To in- 
vestigate this effect, we consider the same synthetic 
GRB as used in H3.2I to model time bins 'a' and 'b' 
of GRB 080916C. We focus on the onset of the GeV 
component, which occurs at i bs,onsct = iobs,tri g + 
0.67 (1 + z) s, where i bs,trig corresponds to the ob- 
server time of the first MeV photons. The MeV pho- 
tons observed at i bs, onset are emitted at radius i?McV 
and it is assumed that the emerging GeV photons 
observed at the same time were emitted by mate- 
rial moving with Lorentz factor TgcV (velocity /3cevc) 
at radius RoeV and time ^GeV (source frame) with 
i Go v ^ Rgcv/c = i bs,onsct/(l + z). The flash of GeV 
photons emitted at i?GeV is assumed to have a power- 
law spectrum with photon slope (3 — —2.2. We define 
a latitude-averaged 77 opacity for GeV photons of en- 
ergy Ebcv by 

e -T 7 ,(£ GoV ) = /e-^^v e °)P(9 c ) 1 ^sine e de e 

/p(e c ) 1 -' 3 sine c de c 

(2) 

where r 77 (i^GoVj ©e) is the opacity seen pho- 
tons emitted at colatitude O c and V (9 C ) = 
(rGeV (1 — /^GcV c °s @e)) is the corresponding 
Doppler factor. The contribution of each colatitude 
to the mean value is weighted by the corresponding 
fluence, leading to the 1 — /3 exponent. We plot 
in Fig. [3] the evolution of the latitude averaged 77 
opacity r 77 at 16 GeV (source frame) as a function 
of i?GcV f° r -RgcV > -RmcV and for different values 
of the Lorentz factor TgcV- When GeV and MeV 
photons are emitted at the same location, we find 
that r 77 < 1 for r Gc v > r G cV,min, same zone — 340, i.e. 
the same limite as in ^3.21 When i?GcV increases, the 
opacity r 77 decreases as expected, which loosen the 
constraint on the minimum Lorentz factor r m j nj GeV 
of the material emitting GeV photons: 



-RGcv/-RMcV 


1 1.2 5.1 13 


Tmin, Go v/^min.GcV, same zone 


1 0.59 0.29 0.15 



This follows approximatively the dependency on 
Re/Ro found in Q (section 2.2.3), i.e. t 77 oc 
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Figure 3: Two emitting region scenario. The 77 

opacity r 77 seen by 16 GeV photons (source frame) 
observed at t obs , onsct — fobs, trig — 0.67 (1 + 2) s (see text) 
is plotted as a function of their emission radius RgcV for 
different values of the Lorentz factor TgcV = 50, 100, 
200, 400 and 800 from top to bottom, [figure from 0]] 



(R e / Ro) 21 " 13 when R c ^> Rq, leading to r m i n .GcV °c 
(R e / Ro) -1 . As shown in §3.21 the detailed modeling 
of the 77 opacity in a scenario where GeV and MeV 
photons are emitted in the same regions leads to a 
reduction of the minimum Lorentz factor by a factor 
~ 2-3 compared to single zone models. The calcula- 
tion presented here shows in addition that the mini- 
mum Lorentz factor can be reduced further more by 
another factor ~ 2 — 8 for FqcV if GeV emission be- 
comes efficient at a radius larger than for MeV pho- 
tons. Assuming that the radiated energy at i?GoV 
is not larger than the radiated energy at RmcV, we 
have checked that the outflow remains optically thin 
for the Thomson opacity due to primary electrons and 
secondary leptons at i?GeV m the case shown in Fig. [3l 
Note that this result on the loosening of the constraint 
on the minimum Lorentz factor does not apply to 
models where GeV photons are entirely due to the 
external shock. Indeed, the small value of t bs, onset 
implies an early deceleration. As the isotropic equiv- 
alent energy of GRB 080916C is huge, this leads to 
a minimum Lorentz factor T > 10 3 in the outflow, 
which is more constraining that the 77 opacity limit. 
The discussion of the effect of a distinct GeV emis- 
sion region presented here is quite simplified and some 
limitations should be kept in mind. If an additional 
GeV component could be firmly identified in GRB 
080916C, the maximum energy i?MeV,max of photons 
associated with the main component should be taken 
into account to derive a new constraint r m ; n on the 



Lorentz factor of the outflow during the MeV emis- 
sion phase. We have assumed here V = 340 as derived 
in H3.2\ using £McV.max = 16 GeV (source frame) but 



r m in will be reduced if E 



McV,max 



is lower. There is 



one further complication: the component produced at 
i?GoV extends probably in the soft gamma-ray range, 
as suggested by the observation of a soft excess corre- 
lated with the high energy component in some GRBs 
such as GRB 090926 [|, GRB 090926 [H and GRB 
090510 It has been assumed here that the anni- 
hilation rate of GeV photons with the seed photons 
produced at RgcV is negligible compared to the an- 
nihilation rate with MeV photons produced earlier. 
This is however not necessarily the case depending on 
the relative intensity of the two components. Clearly, 
a detailed modeling of the emitted spectrum is nec- 
essary to investigate such effects. This is beyond the 
scope of these proceedings and we leave to a forthcom- 
ing study the coupling of the formalism presented here 
to compute the 77 opacity with a detailed radiative 
model such as developed by [l9j . 

5. Conclusions 

5.1 . A new formula of T min 

This study clearly illustrates the need for a detailed 
modeling to constrain the Lorentz factor in GRB out- 
flows. However, when it is not possible, a reasonably 
accurate estimate of r m j n can be obtained from the 
following formula 



[C 1 2 1+2 ^X( ( 8)] ^« 



RmcV 

2 



-Rgc 



1/2 



(! + *)-*= 



(TOeC 2 ) 5 



(3) 

where C\ ~ 4 • 10~ 2 , Ai var is the observed variabil- 
ity timescale, Rqcv/RmcV is the ratio of the radii 
where the GeV and MeV components are emitted, 
and where the high energy spectrum (over a dura- 
tion ~ Ai var ) is assumed to follow a power-law with 
photon index j3 above an observed characteristic en- 
ergy E c : F(E) = F{E C ){E/E c f (ph.cm-^keV" 1 ). 
Energy E max is the observed energy of the most en- 
ergetic detected photons. As usually the spectrum is 
measured over a time interval Ai spcc which is larger 
than the variability timescale Ai var , the normaliza- 
tion F(E C ) entering in Eq.© (fluence at energy E c 
in ph.cm _2 .keV _1 ) must be corrected by a factor 
F(E C ) — F(E C ) x (At var /Ai; S p OC ). This equation can 
be directly applied to Fermi-LAT observations and 
generalizes the usual formula given by 4] by intro- 
ducing two corrections: (1) a more accurate normal- 
ization including a numerical factor C% obtained from 
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the comparison with numerical simulations presented 
in i )3.2l and [7|; (2) the possibility to take into ac- 
count two different emitting regions for MeV and 
GeV photons. The standard limit is obtained with 
i?GoV / ' RyicV — 1 (same region) : then the denominator 
in Eq.(j3|) equals 1. The radius RmcV is estimated from 
the variability timescale by i?MoV — r 2 cAi var / (1 + z), 
which is valid for most models of the prompt emission. 
The radius i?GeV is difficult to constrain without a de- 
tailed model of the high-energy emission mechanism. 
If GeV photons have an internal origin, an upper limit 
for i?GeV is given by the deceleration radius. In the 
future, a measurement of the variability timescale in 
the GeV lightcurve could provide a better estimate of 
this radius. 

5.2. Other effects 

The detailed 77 opacity calculation model pre- 
sented in these proceedings is appropriate and accu- 
rate to study many aspects and consequences of 77 
annihilation in GRBs. In the present work we focus 
on the internal shock model and consider the conse- 
quences and signatures that 77 opacity could have 
in GRB observations, showing that: (i) the tempo- 
ral evolution of r 77 during a burst could favor a delay 
between the MeV and GeV light curves (ii) the 77 cut- 
off transition can be characterized in time-integrated 
spectra. It is usually closer to a power-law steepening 
than to a sharp exponential cutoff. The exact shape 
of the transition strongly depends on the details of the 
GRB dynamics, (iii) for complex GRBs, the 77 opac- 
ity could suppress the shortest time-scale features in 
high energy light curves (above 100 MeV). Only the 
point (i) is discussed in these proceedings: we refer the 
reader to the corresponding paper [7j for more details 
on these different aspects. 
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